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ABSTRACT 

o 

We report spectroscopic observations with the Gemini South Telescope of 

c^ . 5 faint (V~ 20) RR Lyrae stars associated with the Pisces overdensity. At a 

hehocentric and galactocentric distance of ~ 80 kpc, this is the most distant 

substructure in the Galactic halo known to date. We combined our observations 

with literature data and confirmed that the substructure is composed of two 

<^ . different kinematic groups. The main group contains 8 stars and has (Vgsr) = 50 

O ■ km s~^, while the second group contains four stars at a velocity of (Vgsr) = 

^ . —52 km s , where Vgsr is the radial velocity in the galactocentric standard of 

rest. The metallicity distribution of RR Lyrae stars in the Pisces overdensity 

is centered on [Fe/H] = —1.5 dex and has a width of 0.3 dex. The new data 



c/2 , allowed us to establish that both groups are spatially extended making it very 

unlikely that they are bound systems, and are more likely to be debris of a tidally 
disrupted galaxy or galaxies. Due to small sky coverage, it is still unclear whether 
^ ! these groups have the same or different progenitors. 

vn 

o 



Subject headings: stars: variables: other — Galaxy: halo — Galaxy: kinematics 
— Galaxy: structure 



1. Introduction 



X 

^H ' In the last decade, abundant observational evidence has been accumulated in favor of a 



hierarchical formation scenario of the halo of the Milky Way (see a review by lHelmill2008l ). 
The substructures in the halo have been observed both as overdensities in space and as mov- 
ing groups. The tidal streams of the disrupting Sagittarius dwarf spheroidal (dSph) galaxy 
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are the best example of such substructures, with the streams wrapping around most of the 
skv( Ivezic et aliboooi IVannv et allboool : Ivivas et al.lboOll : llbata et alJl2002[lMaiewski et aL 



2003). Ot 



ler kn own substructures include the Virgo Stellar Stream (JDuffau et al. 



Prior et al 



2006 



20091 



1) and sev e ral other small s t ructures associate d with the Virgo overden- 



sitv (INewberg et al.l 120071: IJuric et al. 



(llbata et al. 



2003 



Yannv et al. 



density (IRocha-Pinto et al.l 



2003 



200- 



Orphan stream (IBelokurov et al. 



2008 



Viyas et al.l |2008|), the Monoceros structure 



Ivezic et al.ll2008al ). the Triangulum- Andromeda over 



), the Hercules- Aquila cloud (JBelokuroy et al.ll2007al). the 



2007bl ). and other smaller oyerdensities (IClewlev fc Kinman 



c et al. 



20061 : IStarkenburg et al.ll2009l ) an d streams f Grillmairll2009l : ISchlaufman et al.ll2009[ ). Using 



20001) mult i-epoch data in the stripe 82 region 



Sloan Digital Sky Survey (SDSS; lYor 

(308° < R.A. < 60°, \Dec\ < 1.25°), ISesar et all fl2007D found an overdensity of RR Lyrae 

stars in the halo at a heliocentric and galactocen tric distance of ^ 80 kpc and R.A. ~ 354° 



(termed "Structure J", see Figures 10 and 11 fro mlSesar et al.ll2010l). This structure was sub- 
sequently renamed the "Pisces overdensity" by IWatkins et al.l (120091 ) after the constellation 
where it is located. 

Convincing evidence that the Pisces overdensity is indee d a coherent stru c ture, and not 
a random superposition of stars, was recently presented by iKoUmeier et al.l ( l2009l ). They 
measured radial velocities of 8 RR Lyrae stars located in the Pisces overdensity and found 
that 5 of them have a narrow velocity distribution centered on ~ —75 km s~ in the heliocen- 
tric rest frame. The authors found hints of a second velocity peak in their data, containing 
two or maybe three sta rs, at a lower velocity (< —150 km s~^ in the heliocentric rest frame). 
KoUmeier et al.l (|2009l ) interpreted the main peak as a dwarf galaxy or a disrupting dwarf 
galaxy, and did not consider the second peak as a real group due to a small number of stars. 
In this paper, we expand on their work by reporting radial velocities and metallicities for 5 



RR Ly rae stars in the Pisces overdensity, four of which were not observed by lKollmeier et al. 



( I2OO9I ). By combining these two data sets, we are now in position to confirm the existence of 
a second peak in the velocity distribution, and are able to better explore the spatial extent 
and metallicity distributions of these two substructures present in the Pisces overdensity. 



2. Observations 



2.1. Target Selection 



In Figure [T] we show the spatial distribution of a6-type RR Lyrae (RRab) sta rs in the 
Pisces region. The star s were selected from the publicly available ISesar et al.l (I2OIOI ) catalog. 
KoUmeier et al.l ( l2009l ) obtained spectra of 8 of these stars located in the densest part of 
this region (solid circles and one starred symbol in Figured]). Our targets (starred symbols) 
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for spectroscopic observations were mainly chosen to expand the coverage in direction of 
equatorial right ascension. As seen in Figure [H with the addition of our targets all RR 
Lyrae stars in the densest part of the overdensity now have spectroscopic measurements. 

The positions, periods, and epochs of maximum brightness (HJDq) of RR Lyrae stars 
observed in this work are listed in Tableland their g-hand light curves are shown in Figure [2l 
The periods and epochs of maximum brightness are precise to within a few seconds and 
minutes, respectiv ely, allowing the phase of pulsation to be determined with a 1% uncertainty 
( ISesar et al.ll2010l ). The precise estimate of phase is important when subtracting the velocity 
due to pulsations from the measured radial velocity, as we explain in Se ction [2^ Star 343892 
(indicated by an arrow in Figured]) is in common with lKoUmeier et al.r s sample and provides 
a good check on the consistency of both works. 

The heliocentric distances of RR Lyrae stars, d, are estimated from meas urements of 



their fiux-averaged and dereddened V^-band magnitudes, (V) (for details see ISesar et al. 
20101), and by assuming an aver age value for the ir absolute magnitudes of My = 0.6. This 
value was obtained by using the IChaboyeii (Il999l ) My — [Fe/H] relation 



My = (0.23 ± 0.04) [Fe/H] + (0.93 ± 0.12) (1) 

an d adopting \Fe/H] = —1.5 as the metallicity of RRab stars (the median halo metallicity. 



see 



Ivezic et al.ll2008al ). This metallicity is used when calculating distances even though for 



some stars we measure spectroscopic metallicities later in the paper. The fracti onal error in 



the he liocentric distance due to halo metallicity dispersion {(y[Fe/H] = 0.3 dex, llvezic et al. 
2008aJ) and RR Lyrae evolution off the zero-age horizontal branch, is estimated at 6% (~ 
5 kpc error in distance at 80 kpc). Due to their angular position, the heliocentric and 
galactocentric distances of these stars are virtually the same. The apparent (V^) magnitudes 
of the observed stars range from 20.05 mag to 20.21 mag, corresponding to distances of 78 
kpc and 84 kpc, respectively. The spread in distances is consistent with the uncertainty in 
distance. 



2.2. Spectroscopic Observations 

The spectroscopic observations we re obtained during 2009 Au g 19-21 using the Gem- 
ini Multi-object Spectrograph (GMOS; lAUington-Smith et al.l 120021 ) mounted on the 8.1-m 
Gemini South telescope at Cerro Pachon, Chile. Bad weather on two of the three nights did 
not allow us to observe more than 5 targets. A 1" long-slit was used with the B1200 grating, 
resulting in a resolution oi R = 1820, or ~ 2.5 A at 4500 A. This is equivalent to a resolution 
of ~ 150 km s~^ at H/3. The spectral range extends from 3800 A to 5250 A. The Pisces 
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overdensity targets were observed in two sets of 2700s-exposures centered on Acentrai = 4500 
A and Acentrai = 4550 A (to account for the gaps in the detector), with each exposure set 
bracketed by observations of a CuAr arc lamp. The difference in resolution between the two 
setups (different Acentrai) is very small (2.45 A vs 2.47 A). Nevertheless, each set of spectra 
was reduced and calibrated independently, and only then they were combined. The total 
exposure time (1.5 hours) is < 12% of the pulsation cycle of these RR Lyrae stars. There- 
fore, we do not expect excessive broadening of spectral features due to the changing radial 
velocity during the pulsation cycle (broadening < 13 km s~^). 

Images were processed usin g the qmog packa ge in IRAF. Cosmic rays were removed 
using the L.A. Cosmic routine! ( van Dokumll200ll ). Wavelength calibration was done using 



~ 40 CuAr arc spectral features over the whole spectrum. The median root-mean-square 
(rms) scatter of the fit of the wavelength calibration was 0.03 A, or 2 km s^^. The signal- 
to-noise ratio (S/N) of the final spectra ranged from 20 to 30 at ~ 4000 A. Figure [3] shows 
the spectra of the 5 targets, ordered from the lowest to highest S/N. Balmer lines (starting 
at H/3) and Ca II H and K lines are easily seen in the spectra. 

We observed stars between phases 0.2 and 0.7 of the pulsation cycle in order to avoid 
the discontinuity in the radial velocity curve near the maximum light. As expected, the best 
S/N was achieved for stars observed during the descending part of the light curve (phases 0.2 
to 0.5) rather than during minimum light (phase > 0.5, e.g. star 490555). Radial velocities, 
Vr, were obtained by c ross-correlating our sp ectra with 13 standard star spectra taken from 
the ELODIE catalogs ( jMoultaka et al.l 120041 ) . which were degraded to the resolution of our 



observations (convolved with the instrumental line profile and repixelized) . For the instru- 
mental line profile, we use the line profile measured around the central part of the spectrum 
(i.e., the line profile is assumed to be independent of wavelength). The velocity measure- 
ments are not significantly affected by this assumption, as measured velocities of three HD 
radial velocity standard stars indicate below. The ELODIE radial velocity standard stars 
have spectral types between F2 and F9. The systemic (center-of-mass) v elocity of RR Lyrae 



stars, V^, was obtained following the procedure described in detail in IVivas et al.l ( 12005 



20081 ). Here, we only summarize the most important points. We used the radial velocity 
curve of RRab star X Arietis as a template. The template was shifted in velocity to match 
the measurement for each star at the corresponding phase and the systemic velocity is that 
of the shifted template at phase 0.5. 

Table|2]shows different errors that contribute to the final (systemic) velocity error. First, 



^available at [HREF]http://www. astro. yale.edu/dokkum/lacosmic/ 
^REF]http://atlas.obs-hp.fr/clodic/ 
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there is the cross- correlation error reported by fxcor. This error includes the uncertainties 
due to broadening of spectral features during an exposure, random errors in wavelength 
calibration, and errors due to ELODIE template mismatch. The Ccc value listed in Table [2] 
is the average cross-correlation error value obtained from 13 cross-correlations. The radial 
velocity (K) listed in Tabled] is the mean of 13 Vr values obtained from cross-correlations, 
where each Vr value is weighted by its cross-correlation error. Therefore, radial velocities from 
better template fits contribute more towards the final radial velocity value. The consistency 
of radial velocity measurement obtained from cross-correlations is good (~ 3 km s^^), as 
indicated by the standard deviation crtompiates- 

The cross-correlation error may not fully represent the true error in the radial velocity 
since it does not include systematic errors in the wavelength calibration. To assess this issue, 
we cross-correlated 7 spectra of three radial velocity standards (HD 154417, HD 155967 and 
HD 180482) which were observed during our observing run with the same instrumental setup, 
with ELODIE standards. The velocities of standards were recovered with a mean difference 
of —1.7 km s~^ from literature values; the standard deviation of the differences is only 3 
km s~^. The maximum difference is 6 km s~^, which we assume as the maximum error in 
the wavelength calibration. The error in the radial velocity (cXr) was obtained by adding 
this value in quadrature with the cross-correlation error. The average error of our radial 
velocity measurements is 11 km s~ . Finally, the error calculation for the systemic velocity 
includes ter r as tha t take into account errors in the model radial velocity curve. Following 



Vivas et al.l (120051 ) . we use the following expresion in which the second term accounts for 
likely uncertainties in the phase where the velocity curve passes through the systemic velocity, 
and the third term is related to possible variations in the slope of the velocity curve: 

a^ = al + (119.5 * 0.1)^ + (23.9 * A0)2 (2) 

where A</) = (^ohs — 0.5, and <^ohs is the phase of observation. The contributions of model 
uncertainties, (Jmodeh to final velocities are indicated in Table [2] (cr^ojei — (119.5 * 0.1)^ + 
(23.9* A0)^). These are of the same order as the observed radial velocity errors. Therefore, to 
reduce the uncertainty in the systemic velocity we need to reduce the uncertainty introduced 
by the model, that is, we need to observe the radial velocity curve more than once. The final 
errors in systemic velocities, a^^ range from 10 to 18 km s^ . 

In the case of star 343892 we fi tted the template to both our measurement and the one 



reported by iKollmeier et al.l (120091 ). The fit is shown in Figure HI The radial velocity curve 
of X Arietis fits well to both observations. This good agreement gives us confidence that 
there are no lar g e syst ematic differences between our measurements and those reported by 



KoUmeier et al.l (120091 ). The systemic velocity of this star reported in Table [T] corresponds 



to the one obtained by the fit of the template and differs somewhat from the value assumed 
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by iKollmeier et al.l (see Section [2l3|) . In this case, the final error was calculated as [a, 



a, 



-y, This work) 



-1/2 



(standard deviation of a weighted mean) . 



'y,Koll 



+ 



Finally, we calculated the velocity seen by an observer at the Sun who is at rest with 
r , by assuming solar motion ( 
and vlsr = 220 km s"^ Jsinnev fc MerrifieldlEoOsh . 



the Galactic center, Vg^r , by assuming solar motion of {vu,vv,vw) = (10.0,5.2,7.2) km s ^ 



The spect roscopic metallicities were measured following the method and calibration of 



LaydenI (119941 ) which involves plotting the pseudo-equivalent width of Ca II K line, W(K), 
corrected for interstellar Ca abs orption, against the me an pseudo-equivalent widths of /3, 7, 
and 6 Balmer lines, W(H) (see IVivas et al.l l2005l . l2008l for details). With this method, we 
estimated metallicity errors of 0.15 dex. Due to low S/N, the metallicity for star 377927 has 
higher uncertainty (0.2 dex), and for star 490555 the metallicity was not measured at all (see 
Table [I]). 



2.3. Data from IKollmeier et al.1 (l2009h 



As show n in Figure HI th e re is a consistency between our radial velocities and the ones 
measured by IKollmeier et al.l ( 120091 ). However, the methods that were used to obtain the 
systemic velocity differ between the two works. In order to avoid differences when combining 
samples, we decided to take their reported radial velocities (uncorrected for pulsations) and 
use our method (that is, fitting the template of X Arietis) to uniformly derive systemic 
velocities. When our method is applied to their measurements, we find a systematic difference 
of ~ 15 km s~^(our values are larger than KoUmeier et al.'s). Our r nain conclusions a re not 
affected by this difference. For reference, the original and revised IKollmeier et al.l (120091 ) 
velocities are listed in Table |3l In the following section we use only the revised systemic 
velocities. 



KoUmeier et al. also provide spectroscopic metallicities calibrated using the lGratton et al 



(I2OO4J ) method. For star 343892 that is common to both samples, we measure the metallicity 
of -1.1 dex while KoUmeier et al. measured -1.2 dex. It is quite encouraging that the two 
metallicity estimates are consistent within uncertainties, considering different instruments 
and calibration methods. Given this consistency, we simply combine our and KoUmeier et 
al.'s metallicities when presenting results in the next section. 



3. Results and Discussion 



In this Section we address the following questions: 
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1. What is the mean velocity and velocity dispersion of our sample, and how do these 
values compare to the literature? 

2. Is the observed distribution of velocities a Gaussian or a non-Gaussian distribution? 

3. Are the observed velocity groups bound or unbound systems? 

4. What does the metallicity of stars say about the progenitor(s) of the streams? 

Figure [5] shows the distribution of velocities of 12 RRab stars in the Pisces overdensity 
from the combined data set. The weighted mean velocity, velocity dispersion, and their 
uncertainties for the full sample are (Vgsr) = 7 ± 16 km s~^ and agsr = 56 ± 12 km s~^. If 
the sample is divided into positive (8 stars) and negative velocities (4 stars), the values are 

(YPOS\ ^ 
\ gsr I 

km s""*^. 



50 ± 3 km s-\ alll = 10 ± 3 km s-\ (V;"//) = -52 ± 11 km s~\ and (t"^^ = 23 ± 9 



gsr 



gsr 



The measured velocity di spersion for the full sample is consistent with the halo velocity 
dispersion profile obtained by iBrown et al.l (120 lOl ) using a sample of 910 distant halo stars. 
By extrapolating their Equation 6 to 83 kpc, we obtain a velocity dispersion of 73 km s~^; 
a slightly higher dispersion t han the one we measure, but still within ~ Icr uncertainties. 
Since the (JBrown et al.ll2010l ) sample covers a greater area of the sky than our sample (7300 



deg^ vs. 23 deg^), it is also more likely to contain a greater number of velocity groups at 
large distances than our sample. If these groups have different internal kinematics (mean 
velocities and velocity dispersion), the net effect will be a high er velocity di s persio n in the 
sample with more velocity groups. We have also searched the iBrown et al.l (J2010l ) catalog 
for blue horizontal branch stars that might be associated with the Pisces overdensity and 
found one possible candidate. This star has R.A. =355.575542 deg, Dec=0. 326697 deg, 
Vgsr = 101 ± 17 km s~^, and is at a distance of 74 kpc. The uncertainty in distance is ~ 6%, 
same as for RR Lyrae stars. However, we do not think this star is related to the Pisces 
overdensity as its velocity is at least 3cr away from the positive velocity peak at (Vl^°'^) = 50 
km s~^, where a = a^°^ = 10 ± 3 km s~^. 



gsr 



The question of whether the distribution of velocities shown in Figure is a Gaussian 
or non-Gaussian distribution is an important one, as deviations from normal ity are usually 



interpreted as a signature of velocity group s in the smooth halo component (JHarding et al. 



2001; Duffau et al. 



2006; 



Vivas et al. 



20081 ). T o test for the pr e sence of velocity groups in 



the halo, iHarding et al.l (120011 ) recommend the IShapiro fc Wilkl (119651 ) (SW) statistical test 
of normality to be applied to velocity histograms. This test is sensitive to many different 
deviations from the Gaussian s hape and does not de pend on the choi ce of mean or dispersion 



of th e normal distribution (see iHarding et al.ll200ll Section 5.1 and iD'Agostino &: Stephens 



19861 for more details on the test). 



The null hypothesis of the SW test is that a sample came from a normally distributed 
population. To test a sample, we use the shapiro.test routine provided by the R statistical 
packagcl. This routine evaluates the W test statistic for a sample, and estimates its P-value 
for any n in the range 3 ^ n ^ 5000, where n is the sample size. For example, if the P-value 
of a tested sample is less than 0.05, it can be concluded that the data are not drawn from a 
normally distributed population at the 95% confidence level. 

The P-value is the probability of obtaining a test statistic at least as extreme as the 
one that was actually observed, assuming that the null hypothesis is true. If we know that 
the null hypothesis is true, e.g. we generate samples by drawing values from a Gaussian 
distribution, the test is reliable if the fraction of samples with P < a does not exceed a. 
We use Monte Carlo simulations to estimate the reliability of the P-value provided by the 
SW test for small samples. We generate 10,000 samples and for each sample draw 12 values 
from a Gaussian distribution. Assuming a fiducial value of a = 0.01, we find that less than 
1% of samples have P-values less than 0.01. A similar result is obtained for a = 0.05. These 
results show that the P-value provided by the SW test is reliable even when testing small 
samples: the false rejection rate is less than the adopted value of a. 

The SW test accepts single-valued data points and does not take into account uncer- 
tainties in measurements. To account for uncertainties in measured velocities when applying 
the SW test, we use Monte Carlo simulations. We generate 10,000 samples and for each 
sample draw 12 velocities from 12 Gaussian error distribution (one velocity per Gaussian). 
These Gaussians are centered on measured velocities and have widths equal to uncertainties 
in measured velocities. We then obtain a P-value for each sample using the SW test. The 
most likely (mode) P-value for these samples is ~ 0.01, the average P-value is ~ 0.04, and 
95% of samples have the P-value less than 0.12. This analysis shows that there is only 4% 
probability that the observed sample was drawn from a Gaussian distribution. 

The SW test does not have the ability to assess the rather peculiar nature of the 
shape of the observed distribution (i.e., it cannot distinguish the particular nature of non- 
Gaussianity). To include this information, we consider the small chance (4%) that the 
measured velocities were drawn from a Gaussian distribution, and estimate the probability 
that it would produce the observed distribution seen in Figure O We generate 10,000 sam- 
ples and for each sample draw 12 velocities from a Gaussian centered on 7 km s~^ and 56 km 
s~^ wide. Each sample is then divided into velocities greater and smaller than 7 km s~^, and 
dispersions are calculated for each group {a^"'^ and cr"*^^). In only ~ 3% of generated samples, 
the positive velocities have 7 < 0"^°** < 13 and negative velocities have 14 < a"''^^ < 32, where 



\[nREF]http://www. r-project.org/ 
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the selected ranges reflect the uncertainties in velocity dispersions measured for the observed 
sample. Therefore, even if the observed velocities were drawn from a Gaussian distribution, 
it is unlikely they would create the distribution seen in Figure 13 

In light of these results, we conclude that the hypothesis of normality can be rejected at 
the > 95% confidence level. The next simplest explanation for the observed distribution of 
velocities is that we are observing two velocity groups, one approaching at (V^f^^) = —52 ±11 
km s~^, and the other one receding at (V^^'J.'^) = 50 ± 3 km s^^. 

In Figure [H the stars belonging to two velocity groups are shown in different colors 
(blue for the positive velocity group and yellow for the negative velocity group). One of 
the goals of this work was to probe the extension of these groups in space. Our data show 
that both groups are quite extended in space. The positive velocity group extends from 
349° < R.A. < 357°, while the negative velocity group is mostly centered around R.A ~ 355° 
with one star 10° away, at R.A ~ 344°. At 83 kpc, these angular extensions correspond to 
9 and 15 kpc respectively. The large extensions of these groups are hard to reconcile with 
bound systems. These groups are therefore most likely remnants of disrupted dwarf galaxies, 
such as tidal streams, passing through the stripe 82 plane of obse rvation (similarly to the 



Sagittarius dSph trailing stream, see Figure 12 in lSesar et al.ll2010l for an illustration 



In their work, iKoUmeier et al.l (120091 ) have observed two velocity groups but did not 
consider the negative velocity group as a real system due to a small number of stars (only 
three stars). We have repeated our statistical tests using their sample and original velocities, 
and have found strong evidence of non-Gaussianity in their distribution of velocities. The 
average P-value for their sample is 0.03, and the probability of recreating a distribution of 
velocities with two peaks from a single Gaussian distribution i s less than 0.03%. Therefore, 



we believe there was ample evidence in iKollmeier et al.l (120091 ) for detection of two velocity 



groups, even before four new observations from this work are added. However, the real 
benefit of new observations, and of this work, is the increased angular cover age of the Pisces 



region . The angular extension of velocity groups is now more evident than in lKoUmeier et al. 



( 120091 ). making the velocity groups more likely to be unbound (tidal streams) than bound 



(globular clusters or dwarf galaxies) . 

By combining our and KoUmeier et al. spectroscopic metallicities, we find that the 
positive velocity group has a median metaUicity of [Fe/H] = —1.4 dex and a dispersion of 0.2 
dex. The metaUicity values for this group range from —1.1 dex to —1.6 dex. The metaUicity 
of star 377927 was not used in this estimate due to low S/N of its observed spectrum. 
The negative velocity group seems to be more metal-poor, with a median metaUicity of 
[Fe/H] = —1.7 dex, dispersion of 0.3 dex, and with metallicities ranging from —1.6 dex 
to —1.9 dex. The number of stars in each group is small enough, and the errors are large 
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enough, that we cannot discard the hypothesis that the metaUicities are drawn from the same 
parent population. However, both groups are metal-poor suggesting that the progenitors of 
these two substructures are systems with old, metal-poor populations. 



4. Conclusions 

We confirm the existence of two kinematic groups in the Pisces overdensity, located at 
83 kpc from the Sun, by combinin g new spectroscopic observations with those previously 
obtained by iKoUmeier et al.l ( 120091 ). The spatial extent of groups in the right ascension 
direction suggests that these are not bound systems but rather debris of disrupted dwarf 
galaxies, most likely tidal streams. If these are tidal streams, the question remains whether 
they have the same or different progenitors. The first case (same progenitor) is not an 
improbable one. The Sagittarius dSph galaxy has already shown us examples w here debris 



stripp ed in different epochs overlap in the same region of the sky (see Figure 1 in iLaw et al. 
20051 ) . In such regions large gradients in Vgsr are evident along the same line of sight as some 



streams recede while others approach. The second case (different progenitors) is supported 
by state-of-the-art simulations of galaxy formation. These simulations predict numerous 
substruct ures of different morphologies (streams, shells, clouds) in halos of Milky Way-size 
galaxies ( iJohnston et al.ll2008[ ). Such substructures may overlap and may be difficult to 
separate without wide and deep sky coverage. 



Ongoi ng and upcoming w ide-area surveys, such as Pan-STARRS (JKaiser et al.l |2002[ ) 
and LSST ( llvezic et al.ll2008bl ). will provide the depth and sky coverage needed to separate 
different halo substructures in coordinate space. With single-epoch and co-added observa- 
tions, these surveys will be able to observe main sequence F dwarfs up to distances of 100 
kpc and beyond, and will trace halo substructures with greater fidelity. In addition, detec- 
tions of overdensities in main sequence and RR Lyrae st ars will allow robu st measurements 
of metaUicities using a photometric method described in ISesar et al.l ( l2010l ). 
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Table 1. RR Lyrae Targets 



ID" 



R.A.'' 
(deg) 



Dec*" 
(deg) 



Period 
(d) 



(kpc) 



hjd;^ 



(+2450000 d) 



^■^ ^spectrum 
(+1450000 d) 



Phase-^ 



(km s-l) 



yh 



(km s-i) (km s-i) 



* gsr 

(km s-i) 



[Fc/H]J 



343892 352.469872 -1.171239 0.59731 

377927 350.637200 -1.045226 0.48092 

490555 349.544995 0.856966 0.63091 

719918 344.484556 0.333060 0.62337 

3988771 357.351934 0.759800 0.60580 



83 


3352.59737 


5064.64036 


0.239 


-91 


-57 


10 


58 


-1.1 


83 


3637.73839 


5064.87133 


0.514 


-53 


-55 


18 


63 


-1.7* 


83 


3666.70100 


5065.86722 


0.689 


-49 


-72 


18 


53 




84 


4388.77431 


5064.79280 


0.466 


-186 


-182 


16 


-51 


-1.9 


78 


4412.73725 


5064.72046 


0.226 


-69 


-36 


16 


76 


-1.7 



^RR Lyrac ID number from the lSesar et all 1 12010^ catalog 
'^Equatorial J2000.0 right ascension and declination 
'^ Heliocentric/galactocentric distance 
"^Heliocentric Julian Date of maximum brightness 
"Heliocentric Julian Date when the spectrum was taken 
Pulsation phase when the spectrum was taken 
SRadial velocity (not corrected for pulsations) 

''Systemic (center-of-mass) velocity (corrected for pulsations) and its uncertainty 
' Galactoccntric rest- frame velocity 
J Spectroscopic metallicity, * denotes uncertain values 



hf^ 
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Table 2. Contributions to the error in V^ 



cc ^ templates r model 

(km s^"*^) (km s~^) (km s^^) (km s^^) 



343892 


7.4 


2.3 


9.5 


13.5 


377927 


11.8 


3.7 


13.2 


12.0 


490555 


11.9 


3.2 


13.3 


12.8 


719918 


8.9 


0.5 


10.7 


12.0 


3988771 


6.4 


2.9 


8.7 


13.6 



^Average cross-correlation error 

'^Standard deviation of Vr values obtained from cross- 
correlations 

'^Uncertainty in measured radial velocity 

'^Uncertainty in model radial velocity 
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Table 3. Revision of velocities from iKoUmeier et al.l (J2009[ ) 



ID 


R.A. 


Dec 


Phase 




Vr 




V7 






V-y 




Difference 




Vgsr 










(Observed) 


(Kollmeier 


et al.) 


(Th 


is work) 








(Th 


is work) 




(dcg) 


(dcg) 




(k 


Tl S 


') 


(km s~ 


') 


(k 


tn s~ 


') 


(k 


(11 s 


') 


(k 


ms-i) 


343892 


352.469872 


-1.171239 


0.422 




-63 




-68 






-57 






11 






58 


3115111 


354.116700 


-0.384237 


0.362 




-156 




-156 






-140 






16 






-26 


349151 


354.878980 


-0.157729 


0.436 




-192 




-198 






-184 






14 






-71 


359035 


354.955497 


-0.276307 


0.335 




-192 




-189 






-172 






17 






-59 


3981578 


355.577586 


-0.008393 


0.388 




-82 




-84 






-69 






15 






44 


3974293 


355.600850 


-0.623494 


0.404 




-84 




-87 






-72 






15 






39 


3944324 


355.750750 


-0.173178 


0.425 




-65 




-70 






-56 






14 






56 


3879827 


356.294558 


-0.804956 


0.315 




-73 




-68 






-51 






17 






58 
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Fig. 1. — A zoom-in of Fig. 11 from ISesar et al.l (120101 ) showing the spatial distribution 
of RRab stars located in the vicinity of the Pisces overdensity (red and yellow region). 
Due to their angular position, the heliocentric and galactocentric distances of these stars are 
virtually the same. The colored background shows the observed number density of RR Lyrae 



stars (p) relative to expected smooth model number density {pjrwdei^ ^^- 16 in ISesar et al. 



20101). The observed number density in the red re gion is 10 t i mes g reater than the model 



prediction. The RRab stars were selected from the ISesar et al.l ( 20101) stripe 82 catalog and 



have \Dec\ < 1.25°. The stars spectroscopically observed by iKoUmeier et al.l (120091 ) are 
shown as solid circles, those observed in this work are shown as starred symbols, and stars 
not yet spectroscopically observed are shown as triangles. The arrow points to a star that 
was observed in this work and by iKoUmeier et al.l ( l2009l ). The blue and yellow symbols 
show stars associated with (Vgsr) = 50 km s~^ and (Vgsr) = —52 km s~^ velocity groups, 
respectively. The positions, velocities, and metallicities of spectroscopically observed stars 
are listed in Table [H 
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Fig. 2. — The g-hand light curves of observed stripe 82 RR Lyrae stars, ordere d as in Table [H 
The m agnitudes are corrected for ithe ISM extinction using the map from ISchlegel et al. 
(119981). The light curve s are offset for clarity. The solid lines show best-fit gf-band templates 
from lSesar et al.l (120101 ). 
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Fig. 3. — GMOS-S spectra of the 5 RR Lyrae stars observed in the Pisces overdensity, 
ordered from lowest to highest S/N (top to bottom). An ELODIE template spectrum (nor- 
malized to continuum) is shown at the bottom of the panel. The dashed lines show the 
positions of important spectral features. 
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Fig. 4. — Radial velocity observations of star 343892 (indicated by an arrow in Fig. [T]). 
The solid line is the template of the radial velocity curve of RRab star X Arietis shifted in 
velocity to fit the observations. 
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Fig. 5. — Histogram of Vg^r of the combined sample of RR Lyrae stars in the Pisces 
overdensity. The shaded histogram includes only stars observed in this work. The bin 
size is 20 km s~^, and is slightly larger than the largest radial velocity error (18 km s~^). 
The distribution of observed velocities seems to be bimodal, with two velocity peaks centered 
on (V^sr) = —52 km s~^ and (V^sr) = 50 km s~^, respectively. See the text for a discussion 
of statistical significance of this bimodality. 



